Sampling environmental DNA from trees and soil to detect cryptic arboreal mammals

Environmental DNA (eDNA) approaches to monitoring biodiversity in terrestrial environments have largely focused on sampling water bodies, potentially limiting the geographic and taxonomic scope of eDNA investigations. We assessed the performance of two strictly terrestrial eDNA sampling approaches to detect arboreal mammals, a guild with many threatened and poorly studied taxa worldwide, within two central New Jersey (USA) woodlands. We evaluated species detected with metabarcoding using two eDNA collection methods (tree bark vs. soil sampling), and compared the performance of two detection methods (qPCR vs. metabarcoding) within a single species. Our survey, which included 94 sampling events at 21 trees, detected 16 species of mammals, representing over 60% of the diversity expected in the area. More DNA was found for the 8 arboreal versus 8 non-arboreal species detected (mean: 2466 vs. 289 reads/sample). Soil samples revealed a generally similar composition, but a lower diversity, of mammal species. Detection rates for big brown bat were 3.4 × higher for qPCR over metabarcoding, illustrating the enhanced sensitivity of single-species approaches. Our results suggest that sampling eDNA from on and around trees could serve as a useful new monitoring tool for cryptic arboreal mammal communities globally.


Methods
Our methods consisted of collecting soil and roller eDNA samples at 21 focal trees at two sites; performing eDNA metabarcoding on these samples with a mammal-focused primer set; using a community occupancy model to compare detection rates among species and between collection methods; and assessing the relative ability of these methods to characterize the community of likely arboreal and non-arboreal mammals present. We further compared the ability of a targeted qPCR assay to detect eDNA of one species (big brown bat) in samples relative to metabarcoding. Study area and target species. Field sampling occurred within the Rutgers Ecological Preserve (40.52° N, − 74.44° W; ~ 130 ha) and Morristown National Historical Park (40.77° N, − 74.54° W; ~ 590 ha) in New Jersey, USA. Both sites are dominated by mixed-age second growth Oak-Hickory forest, which is common in the region 34 . We built a database of mammal species occurring at the sites from personal records (unpublished data), park websites, the National Park Service's Inventory and Monitoring Program database, a camera trapping study 35 , and iNaturalist 36 (see Table 1). This database formed the basis for reference for any species detected.
Field sampling. We cleaned and disinfected the sampling pole and soil collection spoons prior to each field use, and transported equipment in sterile bags to avoid contamination (Supplementary Table S1). To test for contamination of field gear, we collected a 'field negative' roller sample upon arriving at the field site (Morristown or Rutgers) by going through the same sample handling procedures described below but without field sampling.
Between July and September 2021, we sampled 21 trees: 12 at Morristown and 9 at Rutgers. Trees were selected on their suspected status as diurnal roosts from a parallel radio-telemetry tracking study involving three species: eastern red bat (Lasiurus borealis; n = 17 trees), northern long-eared bat (Myotis septentrionalis; n = 2 trees), and big brown bat (Eptesicus fuscus; n = 2 trees) (BM, unpublished data). Trees were deciduous (mainly Liriodendron, Quercus, Fagus, and Betula) and 17-81 cm in diameter (mean: 51 cm; see Supplementary Table S2). At 19 trees, we collected soil and roller eDNA within 24 h after a roosting bat was suspected present, and then again on the following two consecutive days (48 h, 72 h), for a total of three rounds of sampling per tree. Two trees received only the first round of sampling for logistical reasons. Thus, we collected a total of 59 sets of roller and soil samples at 21 trees.
We collected soil samples beneath the canopy of focal trees by establishing four transects, oriented 90 degrees from each other, and extending from the trunk out to ~ 3 m, near the outer edge of the canopy. The orientation of the four-way transects was rotated ~ 30 degrees during each subsequent sampling event to avoid overlapping previous sampling locations. Along each transect, we collected a ~ 20 g scoop of soil every 1 m (i.e., 4 per transect), and then combined the ~ 320 g of soil for each tree into a 178 × 305 mm sterile plastic bag. Samples were transported in an ice-filled cooler and stored at − 20 °C until processing.
We collected roller samples on each focal tree trunk following Valentin  www.nature.com/scientificreports/ rolled from the base of the trunk to a height of ~ 3 m around the entire circumference of the tree. Then, we removed the roller from the pole, sealed it in a 304 × 114 mm sterile plastic bag with ~ 400 mL of DI water, and shook the bag while massaging the roller for ~ 30 s to help move any eDNA into solution. We filtered the water using a self-desiccating filter assembly (Smith-Root, Inc., Vancouver, Washington, USA) equipped with a 10-μm polycarbonate track-etched (PCTE) filter membrane and a peristaltic pump. We then used sterilized forceps to remove the filter membrane from the assembly and place it into a 1.5 mL tube filled with 100% non-denatured ethanol. Tubes were stored at − 20 °C until extraction.
DNA extraction and sequencing of environmental samples. For the 59 roller samples, we evaporated ethanol from filters using a vacuum centrifuge, extracted DNA using the DNeasy Blood and Tissue kit (Qiagen), and removed impurities using Ampure XP beads (Beckman Coultier, Pasadena, California, USA). For soil, we extracted DNA from a 5 mL portion of the sample using the DNeasy PowerMax Soil Kit (Qiagen) following manufacturer's protocols. Unfortunately, we were unable to procure enough kits to extract all 59 soil samples due to prolonged shipping delays. To overcome this hurdle, we designed a subsampling scheme, based on haphazardly selecting trees from which to extract DNA from soil samples, until the soil extraction kits in our possession were exhausted. This resulted in a total of 35 soil samples analyzed from 14 trees (see Supplementary  Table S2). We confirmed there was no sampling bias resulting from this subsampling scheme by performing a parallel analysis (see Statistical Analysis and Results). For both soil and roller samples, one extraction negative control was constructed for each extraction batch, consisting of PCR-grade water in place of the sample, but otherwise treated the same. We prepared libraries using the Illumina two-step PCR metabarcoding protocol 37 . For each sample, we ran one initial PCR using the MiMammal-U primer set in replicates of 3 to reduce sample drop-out and PCR bias. The MiMammal-U primer set targets an ~ 200 bp region of the 12S locus and was created to specifically amplify mammal DNA from environmental samples 38 . We included at least one PCR negative control in each PCR run. Reaction preparation for the first PCR followed the Illumina protocol 37 and had the following cycling parameters: 98 °C for 3 min, 35 cycles of 95 °C (30 s), 65 °C (30 s), and 72 °C (1 min), followed by final extension at 72 °C for 5 min, and then hold at 4 °C. For the soil samples, we used the same parameters except with an annealing temperature of 67 °C. The primers also contained partial Illumina adapter sequences to act as binding sites for a second PCR run (described below). We then pooled all the replicates for each sample. After pooling, we performed a bead clean-up for each sample to remove primer dimers and ran a second PCR to add the remainder of the Illumina adapter sequences and add sample-specific dual indexes. For the second PCR, reaction preparation and cycle parameters were identical to those in the Illumina protocol 37 . We performed agarose gel electrophoresis on each sample, and all libraries that produced a visible band were then cleaned with a second magnetic bead clean-up and quantified using the Qubit dsDNA High sensitivity DNA quantification assay. Libraries with no visible band on the gel were not analyzed further, with the exception of negatives, of which all were included. We then pooled all samples in equimolar ratios; 5 μl of each negative were also included in the final sequencing pool. We sent libraries to the Genomics Core Facility at Princeton University for sequencing on an Illumina MiSeq v3 2 × 300nt.

Bioinformatics.
We processed sequence data and made taxonomic assignments following OBITools 39 documentation and Leempoel et al. 23 (see Supplementary Table S3). Briefly, we aligned paired reads, and removed reads with join scores < 40, bases with quality scores < 30, and adapter sequences. We then removed reads with ≥ 21 ambiguous bases, dereplicated reads, removed sequences with read counts < 10 or with length < 80 bp, and filtered all reads for PCR or sequencing errors.
We generated a reference database from all vertebrate sequences available from ENSEMBL 40 and used EcoPCR to match taxa to our sequence reads. For unmatched reads, we assigned taxa manually using BLAST 41 and performed a local BLAST search against a supplemental database of 17 sequences from 9 bat species that we generated from tissue samples (see Supplementary Appendix S1 and Table S4). We then removed all resulting sequences with insert lengths ≥ 200 to exclude non-target taxa (e.g., bacteria) 23 . Additionally, we removed sequences from 10 species that were deemed to be likely contaminants from other projects in our lab (Supplementary Table S5) and five that are commonly associated with contamination in other metabarcoding studies 23,28,30 (Bos taurus, Homo sapiens, Mus musculus, Sus scrofa, Gallus gallus). Following assignment and filtering of taxa, we controlled for contamination by setting a detection threshold for each OTU based on the number of reads that appeared in the negatives, removing any counts that appeared in field, extraction, and PCR negatives from the corresponding samples. Finally, we removed any OTU with a read count < 20 as an additional control for contamination and to reduce noise in the data set.

Big brown bat qPCR methods.
To compare quantitative PCR (qPCR) detection methods with metabarcoding, we re-tested the same 94 extracted roller and soil samples described above using a highly-sensitive qPCR assay for big brown bat that targets a 91 bp segment of mtDNA in the COI region (limit of detection: 2 copies per reaction) 17 . Our qPCR methods matched those of Serrao et al. 17 , except that we used a StepOne Plus Real-Time PCR System (Applied Biosystems, Inc., Waltham, Massachusetts, USA) and lacked an internal positive control. qPCR runs included 3 technical replicates each for all samples, a no-template control, and a six-level standard curve of synthetic DNA (10-fold dilution). Following Serrao et al. 17 , we confirmed species identification in all positive samples by enzymatically cleaning the resulting qPCR product (ExoSAP-IT™ PCR Product Cleanup Reagent, Applied Biosystems, Inc.) and performing bidirectional Sanger sequencing using the same primer set as for the qPCR reaction. The resulting sequences were aligned using Geneious Prime (Biomatters, Inc., Aukland, New Zealand) and identified to species using BLAST as described above. www.nature.com/scientificreports/ Statistical analysis. We assessed the performance of soil and roller eDNA sampling for metabarcoding analyses using a Dorazio-Royle (DR) community occupancy model 42 . These hierarchical models can estimate the effects of covariates on the probability of presence (occupancy) and detection given presence for multiple species simultaneously. Information is shared among species by treating species-specific covariate effects as originating from common distributions. In our model, the dependent variable was the detection (1) or nondetection (0) of mammal species at each focal tree during each sampling event. We included sampling method (soil or roller) as a covariate in the detection sub-model to evaluate differences in detection probability between the two methods. In addition, we visualized mammal community data using Venn diagrams, taxonomic heat trees 43 , species accumulation curves, and non-metric multidimensional scaling (NMDS) ordination plots 44 .
We compared the performance of qPCR versus metabarcoding methods for detecting big brown bat eDNA using multimethod occupancy models 45 . These multilevel models use joint detection information (e.g., from qPCR and metabarcoding) to inform the availability for detection, θ, and detection probability given availability, P(detect | θ). In this case, availability refers to the probability of big brown bat eDNA presence within a sample. The dependent variable was the detection (1) or non-detection (0) of big brown bat by each method. We used a covariate for sampling method (soil vs. roller) within the availability sub-model, and a covariate for molecular approach (qPCR vs. metabarcoding) to estimate the method-specific probability of detecting big brown bat eDNA given availability within the sample. We then estimate the total 'per-visit detection rate' of bats at focal trees by multiplying θ by P(detect | θ). For both the DR and the multimethod occupancy model, we used the per-visit detection probability (p) to estimate the cumulative probability of detecting a species at least once given n repeated samples using the formula 1− (1 − p) n (e.g., Ref. 24 ).
All models were fit in a Bayesian framework with non-informative priors (see Ref. 42 , pp. 607 and 690) using JAGS and jagsUI in program R 46,47 . We ran 3 chains of 30,000 iterations each, including a burn-in period of 10,000, keeping every 10th draw. Model convergence was assessed by examining trace plots and Gelman-Rubin statistics (r-hat < 1.1). We compared estimates of per-visit detection probability for the various methods by examining and plotting posterior predictive distributions.

Results
We performed roller and soil eDNA sampling at 21 trees, obtaining three samples per tree at all except for two trees, which each received a single sample. We analyzed all 59 roller samples from 21 trees, and a subset of 35 soil samples from 14 trees (1-3 samples per tree; mean = 2.5/tree). Forty-seven of 59 roller samples and 35 of 35 soil samples had measurable read counts with a total of 9,852,404 and 4,731,187 reads matched to roller and soil samples, respectively. After all filtering steps (see Supplementary Table S6), contaminant reads still remained in 5 of our 46 negative control samples, all associated with two taxa: Canis lupus familiaris (7702 reads of 4 OTUs) and Glaucomys volans (a single read of 1 OTU). After removal of these negative control reads from samples, 2,548,768 total reads across 46 samples remained for roller samples, and 33,314 total reads across 23 samples remained for soil samples (Supplementary Table S6).
A total of 359 OTUs remained after all filtering steps, attributable to 16 mammal species: 8 arboreal species (5 detected with soil, 8 with roller) and 8 non-arboreal species (3 detected with soil, 7 with roller; Figs. 1 and 2). We identified 28 total species (14 arboreal, 14 non-arboreal) as likely to be present at the sites based on our list compiled from external sources in addition to our metabarcoding results (Table 1, Fig. 1A). Two of those species (house mouse and little brown bat) were not available to be detected using metabarcoding methods as they were excluded due to contamination concerns (Table 1, Supplementary Table S5), and one species (American black bear) was not available for detection because the MiMammal-U primer set is not effective for bears 38 . Thus, the 16 species we detected using metabarcoding represents 64% of the 25 total species expected to be present and available for detection. Of the 9 species that were not detected using metabarcoding, 4 were considered arboreal, all of which were bats (Table 1, Fig. 1A). We detected 2 of 6 expected bat species: big brown bat, a tree-roosting species, and eastern small-footed bat, a species that roosts in rock crevices. Four species (brown rat, vole sp., domestic dog, and small-footed bat) were detected by metabarcoding but did not occur in our compiled list of likely species (Table 1, Fig. 1A).
In general, higher read counts were found in roller versus soil samples and for arboreal versus non-arboreal species. For roller samples, the mean number of reads per sample averaged 4849 ± 2094 (SE) across species for arboreal mammals (range = 56-14,699, n = 8), and 537 ± 353 for non-arboreal mammals (range = 0-2901, n = 8; Fig. 1B). For soil samples, the number of reads averaged 82 ± 48 for arboreal mammals (range = 0-349, n = 8), and 40 ± 33 (range = 0-270, n = 8) for non-arboreal mammals. The number of reads for bats averaged 59 ± 3 (range = 56-61, n = 2) for roller sampling and 14 ± 14 (range = 0-28, n = 2) for soil sampling. The discrepancy in mean read depth observed between sampling methods is partly explained by the removal of a large number of microbial reads from soil samples during bioinformatics (> 95% of reads removed for soil compared with ~ 50% of reads for roller; Supplemental Table S6).
Despite having fewer samples and fewer reads per sample, our soil method resulted in a generally similar characterization of mammal community composition to the roller method (Figs. 1 and 2; see also NMDS ordination plot, Supplementary Fig. S1). Species accumulation curves (Fig. 3) revealed that the number of species detected had not yet leveled off with either method, suggesting that sampling additional trees likely would have added more species. In addition to mammals, we also detected seven other vertebrates (six birds and one salamander), all of which were detected using roller methods, but only one of which (American Robin, Turdus migratorius) was detected using soil methods (Fig. 2).
Arboreal mammal species with the highest detection probability according to the community occupancy model were also the species with the highest read counts (Figs. 1B and 4 www.nature.com/scientificreports/ non-arboreal white-tailed deer also showed relatively high detection probability estimates within both roller (0.58) and soil (0.23) samples (Fig. 4). Based on these estimates for southern flying squirrel, the most cryptic non-bat arboreal mammal that we detected, ~ 4 roller sampling visits to an 'occupied' tree would be needed to ensure 95% confidence of detection, compared with ~ 12 visits required with our soil eDNA methods (Fig. 5).  For big brown bat, metabarcoding detections occurred in 3 of 94 samples (all roller samples) at 3 of 21 trees, while qPCR detections occurred in 10 samples (9 roller samples and 1 soil sample) at 6 trees. Together, metabarcoding and qPCR detected big brown bat DNA in 11 samples at 7 trees, including one of the two trees identified as likely containing a roosting individual of this species. The multimethod occupancy model revealed that detection probability was 3.4 × higher using qPCR-based detection methods (49%, 95% CI = [21%, 86%]) than when using metabarcoding methods (15%, [3%, 41%]). Based on these estimates, ~ 7 visits to an occupied tree would be required to establish presence of big brown bats with 95% confidence using qPCR methods, while ~ 27 visits would be required using metabarcoding methods (Fig. 5).

Discussion
Fully terrestrial eDNA sampling approaches offer a potentially powerful addition to biodiversity monitoring efforts 23,24 . However, protocols for using eDNA-based methods to characterize terrestrial biodiversity, and vertebrate communities in particular, are still nascent 27,28,30 . In this study, we show for the first time that an eDNA metabarcoding approach can be used to broadly characterize tree-dwelling mammal communities by sampling tree trunks and surrounding soil. Our findings add to recent work (e.g., for reptiles 22,24 ) showing that surface eDNA collection methods, which are relatively untested compared with soil-based eDNA methods, can also be effective at detecting terrestrial vertebrates. Further, we demonstrate that supplementing metabarcoding detection with qPCR-based methods can greatly improve sensitivity, a potentially important consideration for monitoring schemes focused on rare taxa (e.g., Refs. 11,12 ). Together, our results have significant implications for global biodiversity conservation as the broader guild of arboreal vertebrates includes highly threatened 5,48,49 , as well as invasive alien species 50 , that are often cryptic, inhabit inaccessible locations, and are therefore challenging to monitor. Our methods captured over 60% of the mammalian diversity expected at the sites, and a similar fraction of the subset of arboreal species, despite sampling only 21 trees. Species accumulation curves suggest that more species would likely have been added with increased sampling effort. These results broadly agree with those of Leempoel et al. 23 who found that soil eDNA metabarcoding well characterized mammal communities in California chaparral. However, in both our study and that of Leempoel et al. 23 , some conspicuous absences were evident. Bats comprised all of the arboreal species that we expected but failed to detect at our sites using metabarcoding (Fig. 1A). Leempoel et al. 23 also noted a lack of bat detections (2 of 14 possible taxa detected), which they suggested could be due to low efficiency of either the 12S primer set or of their soil sampling methods for that order. While both reasons could also apply to the lack of bats detected in our study (discussed further below), the performance of the 12S primer set very likely contributed to our lack of American black bear detections as MiMammal-U primers are known to be ineffective at amplifying bear DNA 38 . These challenges highlight the reality that false negatives and varying detectability among species are common issues to all survey approaches, including eDNA metabarcoding. Our study represents a rare example among metabarcoding studies in that it uses repeated sampling and community occupancy models to quantify false negative rates. This quantitative approach, coupled with continued experimentation with different molecular techniques and survey methods (e.g., Refs. 23,27 ), will be vital to helping researchers decide how eDNA metabarcoding methods will fit into existing biodiversity monitoring efforts moving forward. www.nature.com/scientificreports/ Although our results suggest that sampling for tree-roosting bats using eDNA metabarcoding still requires further research and optimization, our approach likely has application to characterizing communities in a much broader range of arboreal species globally. Geographic regions with multiple elusive arboreal mammals of management interest-for example, gliders and tree kangaroos in Australasia, or primates in the global tropics-may be particularly suited for a metabarcoding approach for community-level assessments 4,8,9,49 . It may be especially useful for rapid biodiversity assessments (RBAs 51 ) in remote forested environments, where the ability to collect multiple samples relatively rapidly without regard to time of day would be a key advantage 27 . Existing survey methods to monitor arboreal mammals tend to be optimized for particular groups of species, often segregated by body size and behavior, with no suitable single method available to characterize all members of the guild 4,8,9,16,49 . Diurnal and nocturnal species, for example, often require separate survey methods or timing 8 . While camera traps capture both diurnal and nocturnal species, they typically miss smaller species 16,23 . The need for multiple methods to survey for nocturnal and diurnal, or large and small, species separately raises the cost of sampling and can result in datasets that are difficult to compare across sites because of inherent sampling biases 8 . Excluding bats, we found encouraging results for both diurnal and nocturnal arboreal species of a broad range of body sizes, detecting all seven expected species (Fig. 1A). While more work is needed to assemble robust genetic reference libraries before global arboreal mammal monitoring with eDNA metabarcoding will be broadly feasible, a clear advantage of the technique remains the power to detect a broad swath of species, with widely varying morphologies and behaviors, with a single method 23,27,28,51 .
The promise of eDNA metabarcoding approaches for at least some arboreal guilds is well illustrated by our results for southern flying squirrel, Glaucomys volans. Like other flying squirrels (Tribe: Pteromyini), this species is strictly nocturnal, highly arboreal, and tends to get injured in live traps, making it difficult to directly observe and monitor 48,52 . Yet G. volans eDNA was readily detectable using metabarcoding in our study, occurring in 19-26% of soil samples and 47-52% of roller samples across both sites. Our similarly encouraging results for detecting other squirrels (Sciuridae) also bode well for management applications. For example, the methods would enable fine-scale mapping of habitat use in places such as the United Kingdom where native red squirrels (Sciurus vulgaris) are outcompeted by eastern gray squirrels, or the Delmarva peninsula (USA) to support the conservation efforts for the Delmarva fox squirrel (Sciurus niger cinereus) 53 . Further research is needed to determine the extent to which our results for squirrels generalize to other taxa with similar active tree-climbing lifestyles (e.g., gliders 4 , primates 49 ).
Our finding that soil samples revealed fewer species, had lower detection probability, and had lower read counts than roller samples, even for some non-arboreal species like white-tailed deer, likely reflects multiple factors. First, soil and tree bark represent markedly different biological and chemical environments that likely Figure 5. The cumulative probability of detecting a southern flying squirrel (left plot) or big brown bat (right plot) at least once at an occupied tree given increasing numbers of sampling visits. The plot, for southern flying squirrel, a cryptic arboreal species, illustrates the differences between tree bark surface ('roller') versus soil eDNA collection methods within a metabarcoding detection framework. The plot for big brown bat illustrates the difference between qPCR and metabarcoding molecular detection methods on samples collected using roller eDNA techniques. Shaded areas show 95% credible intervals. The gray dashed horizontal line represents a cumulative 95% certainty of detecting the species at least once. www.nature.com/scientificreports/ differ in eDNA quantity by species, eDNA persistence rates 54,55 , and microorganism abundance. The latter may be especially pertinent to our study as we observed a relatively large drop in the number of reads after removal of microorganism reads, especially for soil samples. This suggests that performing additional purification steps prior to sequencing could boost the ability of both methods, and especially soil eDNA, to detect target species by increasing mammalian sequencing depth. Other in-lab factors, such as method of extraction 23 or choice of primers, similarly have the potential to influence the recovery and amplification of target species' DNA and should be the focus of future research. Next, our focal trees were not chosen to occur near any special attractants or areas of multi-species use, such as saltlicks or water sources, which has proven successful in other vertebrate eDNA studies 18,25,31,32,56 . It is possible that adding a broader range of soil sampling sites, including some targeted towards other guilds (e.g., burrow users 32,56 ), would have yielded a more complete inventory. Nevertheless, both soil and surface methods have advantages over the much more commonly-used metabarcoding approaches that rely on natural water bodies for assessing mammal communities 16,17,[27][28][29][30]38 as they are not limited to where these features occur. Our study is the first to suggest that surface eDNA metabarcoding methods can be a powerful supplement to established soil-based methods of characterizing mammal communities, especially for arboreal species.
As noted, bats were especially lacking from our eDNA metabarcoding results, with only two of six likely species detected. Notably, our metabarcoding species list lacked two of the bat species that our sampling scheme was designed around (eastern red bat and northern long-eared bat) and for which we had confirmed recent presence at the sites ( Table 1). The lack of northern long-eared bat detections may directly relate to recent precipitous population declines (~ 99%) caused by white-nose syndrome 57 . However, the lack of eastern red bat detections was especially surprising as roosting of this species was suspected based on telemetry in 17 of our 21 target trees. Reasons for this omission may relate to the fact that eastern red bats roost singly on small twigs and in leaf clusters, and therefore may not leave much DNA on tree trunks. Another possibility is low efficiency of the 12S primer set for bats, although we were unable to find information about this in the literature. It is notable that Leempoel et al. 23 had a similarly poor representation of bats with comparable soil-based methods. However, our metabarcoding results did indicate that we are capable of detecting even uncommon, or at least unexpected, bat species with our methods. Eastern small-footed bat, which is typically viewed as a rock-roosting species and is considered endangered by the International Union for Conservation of Nature (IUCN) 58 , was detected in both soil and surface eDNA samples from Morristown National Historic Park. This species was not otherwise confirmed as present at the site until a year later, in spring 2022, when it was caught in a mist net (BM, unpublished data). Our results with respect to bat detections, along with those of others 23 , underscore the need for further research to adapt eDNA metabarcoding methods to this vulnerable group, which could contribute much needed demographic and distribution information. This is especially urgent as 18% of bat species are listed as "data deficient" by the IUCN, while 57% lack basic population trend information 5,58 .
Our comparison of qPCR to metabarcoding detection methods for big brown bat represents a hopeful result for the use of eDNA to monitor rare vertebrates that are of particular conservation interest. It is well-known that qPCR-based eDNA surveys targeted towards individual species return higher detection probabilities and have greater power at low abundance, than metabarcoding approaches 59 . Our results agree, showing for the first time that adding a qPCR step in the analysis of surface and soil eDNA samples can be effective for detecting bats in forested environments. The addition of a qPCR step opens the door for developing species-specific assays to increase detection power for endangered or elusive bat species, or other cryptic arboreal mammals 49 . Emerging molecular detection approaches such as droplet digital PCR have the potential to increase this sensitivity even further 59 . Like other eDNA-based tools and survey tools in general, careful consideration of sampling effort, the natural history of target species, and the configuration of different field and molecular methods will be key to optimizing our approach to characterize mammal communities, or to target a particular species, in different regions.
Although eDNA surveys are not inexpensive given the need for both fieldwork and molecular analyses, they can be cheaper than conventional approaches, especially if such approaches require many hours of fieldwork or expensive equipment 60 . Thus, the relative cost-effectiveness of surface or soil eDNA surveys will depend heavily on the mammal communities of interest, the mix of methods that must be employed to effectively sample them, and the purpose of the sampling efforts. However, even if costs are increased, eDNA surveys can reduce field time to the extent that they can improve detection rates, either by replacing or supplementing conventional sampling methods (e.g., as a supplement to visual observations). With higher detection rates, fewer visits are required to achieve the same results. This operational efficiency would be especially advantageous when field conditions present safety risks, are intrusive to sensitive habitats, or are challenging to access. For example, adding surface eDNA sampling to existing visual surveys of eastern wood rats (Neotoma floridana), a cryptic mammal that inhabits steep, rocky slopes in the eastern US, could likely increase detection power, thereby reducing the need for additional risky and costly sampling visits. More studies involving direct comparisons among methods (e.g., Refs. 23,24,30,60 ), in a variety of ecoregions, are needed to determine the extent to which incorporating our methods into existing vertebrate monitoring workflows would increase efficiency.
Finally, we detected other vertebrates, including seven birds and one salamander, in soil and surface eDNA samples, despite our use of a mammal-specific primer set. This is similar to results from Leempoel et al. 23 in California using the same primer set, in which six bird species were detected. We found that surface eDNA detected more bird species than soil, perhaps for the same reasons as for mammals (above). Our results provide evidence that surface eDNA surveys, with taxon-specific primers, could be used to survey bird communities, or used to target particularly rare species in forested ecosystems (e.g., Ref. 61 ). Our detection of a salamander, coupled with recent promising research into reptile detection using surface eDNA methods 22,24 suggests a broader potential for applications with other vertebrates as well. Finally, both surface and soil eDNA metabarcoding can be expanded beyond forests, providing insight into their effectiveness in other habitats (e.g., caves 17 or talus slopes). Our study